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Recently, Fukushima [Phys. Rev. B 78 115105 (2008)] proposed 
a systematic derivation of the GutzwiUer approximation for the t- 
J model. In the present paper, using this approach we construct 
an effective single-particle Hamiltonian, which leads to a renormal- 
ized mean-field theory (RMFT). We also use the method proposed by 
us recently and based on the maximum entropy principle (MaxEnt), 
which in turn, yields a consistent statistical description of the problem. 
On the examples of non-magnetic superconducting d-wave resonating 
valence bond (dRVB) and normal staggered-flux (SF) solutions, we 
compare two selections of the Gutzwiller renormalization schemes, i.e. 
the one proposed by Fukushima with that used earlier by Sigrist et 
al. [Phys. Rev. B 49, 12 058 (1994)]. We also confront the results 
coming from our variational solutions with the self-consistency con- 
ditions build in, with those of the non-variational approach based on 
the Bogoliubov-de Gennes self-consistent equations. Combination of 
the present variational approach with the new renormalization scheme 
(taken from Fukushima's work) provides, for t/J = 3, an upper criti- 
cal hole concentration Xc ~ 0.27 for the disappearance of the d-wave 
superconductivity. Also, the hole concentration x « 0.125 is obtained 
for the optimal doping. These results are in rough accordance with 
experimental results for high-Tc superconducting cuprates. 



PACS: 05.30.-d, 71.10.Fd, 75.10.Jm. 

1 Introduction 

The t-J model [Ij is commonly regarded as a minimal model capable of de- 
scribing correctly the essential physics of the cuprate high-temperature su- 
perconductors. In its simplest form, used in the present paper, it is expressed 
by the following Hamiltonian 
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Here Pg = Ylii^ ^ ^it'^j;) is a Gutzwiller projection operator, ensuring that 
no doubly occupied sites are present. 

Unfortunately, the rigorous treatment of the t-J model is limited to very 
special cases. This is both due to the interaction term (Sj ■ S-,), as well as 
to the presence of Pq operators. Even in the J = limit, Hamiltonian ([1]) 
is not an independent-particle one. To proceed further, ([T]) may be treated 
within the mean-field approximation. This procedure yields an single-particle 
Hamiltonian H, being the mean- field analogue of Ht-j ([I]). The resulting 
effective description is termed the renormalized mean-field theory (RMFT) 

The crucial point in the construction of RMFT is an approximate treat- 
ment of the Gutzwiller projection. In a broader perspective, this problem is 
closely related to an analytic evaluation of the expectation values of operators 
with respect to the following variational state 

m = PgcI^o) = n4"'^rf(l - n,tn4)l^o) (2) 



In the above, Pq was replaced by a more general Gutzwiller correlator Pgc, 
differing from Pq by the presence of the so-called fugacity factors Xi^ [3]. 
For Ajo- = 1 we recover Pg- I^&o) is an uncorrelated single-particle state, 
which, within the framework of RMFT, is chosen as the eigenstate of the H. 
Explicitly, we are interested in evaluation of the expressions like 

O is an arbitrary operator, {Ag} = Ai,A2, . . . ,Am (= ^4) are the relevant 
mean-fields, i.e. the expectation values of the corresponding single-particle 
operators {Ag} = Ai, . . . , Am, and both H and |\E'o) are usually A- depen- 
dent. Consequently, fo{A) is a C- valued function of the mean-fields. Each 
prescription of the form will be termed the renormalization scheme (RS). 
Note, that usually the RS ([3]) may be given a more specific form 

^6)a^'^^fiA)'^^^f(Am. (4) 

In the above, g'^{A) is termed renormalization factor. 

There are many different renormalization schemes ([3]) proposed in the 
literature. The simplest take into account only on the local (site-dependent) 
mean fields (e.g. local charge density or magnetization), [2]-[Z] The more 
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advanced include also the inter-site quantities (mean fields defined on bonds, 
e.g. pairing amplitude), cf. the Refs. [3], [H]-[ID]- In particular, the scheme 
of Ref. [8] has been used in Refs. [TT] and [12] and also in our previous work, 

m- 

In this paper we focus our attention mainly on the renormalization scheme 
of. Ref. [3], which, in our opinion, is the most promising one devised so far. 
First, due to the presence of fugacity factors in ([2]), the Gutzwiller projec- 
tion does not change the local densities of spin up (down) electrons, i.e. 
{nia)G = (nia) = nitj. Next, this formalism may, in principle, be systemati- 
cally extended beyond the second order in the inter-site quantities (formulas 
of Ref. [3] are provided up to this order). Also, it can be relatively easily 
applied to the extensions of Hamiltonian ([T]), or to the more complicated 
symmetry - breaking situations. 

Suppose, that we have constructed the RMFT specifying appropriate 
single-particle mean-field Hamiltonian H. Next, we have to determine the op- 
timal values of mean-fields appearing in a model. One possibility is an appli- 
cation of the non- variational self-consistent approach based on the Bogoliubov- 
de Gennes (BdG) equations. However, this route, apart from its other draw- 
backs, [H], encounters a serious difficulty in the following sense. Namely, the 
renormalization scheme of that Ref. [3] is in general not of the form (jl]). This 
feature does not allow for an unambiguous identification of the renormaliza- 
tion factors and hence the construction of the effective MP Hamiltonian. 
Namely, there is no clear way to ascribe the corresponding single-particle 
operator expression to the expectation values ([3j). Even if some way of re- 
placing (O) ([3]) by the operator counterpart may seem more natural than 
others, there is in fact no unique way of carrying out such procedure. Also, 
different such ways obviously yield different versions of RMFT Hamiltoni- 
ans and consequently different BdG self-consistent equations, hence different 
predictions of the model. 

This lack of uniqueness does not appear within the variational method 
proposed by us recently [Hj, as will be discussed in detail below. 

The main aim of the present paper is twofold. First, we compare different 
renormalization schemes within the RMFT effective single-particle picture for 
the t-J model. Moreover, as in to Ref. [13], we compare two distinct methods 
of solving such MF models: the present variational approach (labeled as var) 
and a non- variational one, based on the Bogoliubov-de Gennes self-consistent 
equations (labeled as s-c). 

The paper is organized as follows. In Section [2] we present those parts of 
our approach which are relevant to the present discussion. In Section [3] we 
present our numerical results, first for the d-wave superconducting (dRVB), 
(13.11 and 13.21) . and then for the staggered-fiux (SF) solutions (13.31) . respec- 
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tively. In Appendix A we provide some of the technical details omitted in 
main text. Section H] contains a concluding remarks. 

2 Formalism: the method 

2.1 Self-consistent variational approach to the mean- 
field models 

The most of standard mean-field models may be solved by employing two 
methods. First, we may invoke the variational procedure, i.e. minimization 
of the MF grand potential or ground state energy, with respect to the values 
of the order parameters (mean fields). Secondly, by using the self-consistency 
(Bogoliubov-de Gennes) equations, expressing the basic fact that the mean- 
fields are averages of the corresponding operators. However, those two routes 
are equivalent only for MF Hamiltonians of the Hartree-Fock type (e.g. the 
BCS Hamiltonian) [H [IS]. In the case of RMFT for the t-J model, the 
self-consistency of the MF formalism may be spoiled by unwary application 
of such variational procedure, for both nonzero temperatures and for T = 0. 
This due to the non-Hartree character of the RMFT Hamiltonian, caused by 
MF treatment of the Gutzwiller projection. 

In such situation, the non variational self-consistent method based on the 
Bogoliubov-de Gennes (BdG) equations is usually applied to those advanced 
versions of RMFT, [H El [71 [HI [12]. However, this route suffers from serious 
drawbacks [H]. 

The other solution is the appropriate modification of the variational 
method. Examples of such self-consistent variational MF approach may be 
found in Refs. [SI [SJ [HI [ID], [IS]- [IB]- The approaches of those References 
exhibit various degree of generality, and differ from each other with respect 
to technical details. Yet, all are based on the variational principle of quan- 
tum mechanics (minimization of the expectation value of the Hamiltonian), 
applied to the MF case. As such, they are applicable to T = situation only. 

Our point of departure is different. We base our approach ([131 [Hj) on 
the maximum entropy principle (MaxEnt) |tl9j, which is the basis of Bayesian 
mathematical statistics. In particular, it may be viewed also as the basis of 
the standard (non mean-field) statistical mechanics [TH [20] . The application 
of MaxEnt inference to the non-standard case of MF statistical-mechanical 
description is a natural extension of this fundamental principle ([14j, c.f. 
also 121J)- To ensure the self consistency of MF description, it is necessary to 
introduce additional constraints, not present in the standard statistical me- 
chanics. This could be achieved in the most natural manner by the Lagrange 
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multiplier method. Explicitly, starting from an arbitrary MF Hamiltonian 
H{A), we define yet another MF Hamilton operator Hx according tc0 

M 

H ^ Hx = H-J2^s{A-A,). (5) 

s=l 

The correct grand-canonical MF density operator is then given by 

PA = exp ( - /3{Hx - /iiV)), Zx = Tr[exp ( - P{Hx - /iiV))]. (6) 

Next, the generalized grand-potential Landau functional is defined as 

T{A,X) = -/3-Mn(Tr[e-^(^^-'^^)]) = -/3-Mn(ZA), (7) 

with the inverse temperature /3 = l/ksT . The equilibrium values oi A = Aq, 
A = Aq are the solution of the set of equations 

Va^ = 0, Va^ = 0, (8) 

for which J^{A, A) has the lowest value. By taking the derivatives with respect 
to A only, and subsequently putting A = 0, 

Va-F = 0, A = 0, (9) 

we obtain standard Bogoliubov-de Genes (BdG) self-consistent equations. In 
such a case we denote the chemical potential /i = ^s-c, which corresponds to 
ft = fi + X for the var method. The thermodynamical grand potential fl and 
the free energy F are defined respectively as 

n{T, V, fi) = J^{T, V, /i; io(T, V, /i), Ao(r, V,fi)), F = Q + fiN (10) 

for var, and similarly (with fig-c, Aq = 0) for s-c method. 

Obviously, the present formalism is valid for any non-zero temperature, 
but not for T = 0. Consequently, we should replace pure states {ipo) in Eqs. 
([2]) - (gD by mixed ones px- 

On the other hand, the Gutzwiller approximation in a form ([3]) or (jl]) is 
devised for a ground state, hence for temperature T = 0. However, the solu- 
tions of the mean field models (and the RMFT of the t-J model in particular), 
obtained for non-zero, but sufficiently low T are practically identical to those 
of the real T = analysis. Consequently, px is practically indistinguishable 

^Lagrange multipliers appearing in ([S]) should not be confused with the fugacity 
factors Xia- appearing in Eqn. ^ 
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from \^po), and thus the apphcation of the finite-temperature formalism in 
the /3 — oo hmit is fully justified. 

Reader interested in details of our method may consult Ref. [H]. Below 
we recall some of its features, which are important from the point of view of 
the present application. They are justified in the Appendix A. 

First, the A - dependence of the nontriviaH part of the Hamiltonian Hx{A) 
^ is completely determined by the set of the single-particle operators out 
of which Hx{A) and H{A) are composed. 

Moreover, two MF Hamiltonians H^^'' and H'^'^^ (or H^^"^ and H^^'^), hav- 
ing the same nontrivial operator part and the same A -dependence of the 
expectation value, i.e. {H^^) = {H^'^^) = {H^^) = (-ff^^^), are equivalent. 
Namely, they yield identical equilibrium values of mean fields, quasi-particle 
energies, thermodynamic potentials etc. This feature guarantees that, when 
constructing RMFT Hamiltonian, different assignments of the operator ex- 
pression to the expectation values f|T3|) lead to identical results. 

Parenthetically, the values of the Lagrange multipliers A may differ be- 
tween two such equivalent MF Hamiltonians This indicates, that 
A have no unambiguous physical interpretation by themselves, but only in 
the certain combinations with the mean fields. Only the Lagrange multipli- 
ers related to the quantities of a priori known average values, e.g. particle 
number A^, are identical for all such equivalent Hamiltonians. 

2.2 Application to the renormalized mean- field t-J model 

For the RS that can be given the form (jlj), with Xi„ = 1 in ([2]), and in the ab- 
sence of the long-range antiferromagnetic order, the mean-field Hamiltonian 
H may be taken in a form [H [Zl El [IH [I2] 

3 

^ = Yl i'^ijdljclcj^ + H.c.) - -JijgfjiXjiclcj^ + H.c. - \x^j\^) 

- Ei'^^^-4(^^.sW-. + H.c.-|A,,n. (11) 
(«i>o- 

In the above expression, (cjo-) are ordinary fermion creation (annihilation) 
operators, Xij = {4aCja), and Aij = {ci_„Cj„) = {cj_^Ci^) are respectively, 
the hopping amplitude (bond-parameter) and the RVB gap parameter, both 

^i.e. not proportional to the unit operator 
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taken for nearest neighbors (ij). (O) denotes the average value of the oper- 
ator O evaluated with the help of a MF statfi (m. 

In (fTTj) . the interaction term (Sj ■ Sj) of Hamiltonian ([1]) has been treated 
within the Hartree - Fock decoupling in both the particle - particle (Ajj) and 
particle - hole (Xij) channels. The renormalization factors gjj and gfj, are 

given by (gD, with 6 = Y.(j^j)„{g\jc\aCj<j + H.c.) and 6 = Si- S^, respectively. 
Their explicit form depends on the approximation used to obtain r.h.s. of 

©• 

As mentioned above, the non-Hartree-Fock character of RMFT Hamilto- 
nians H is due to MF treatment of the Gutzwiller projection. Consequently, 
in order to obtain Hx, we must add to H the constraints corresponding to 
the mean fields appearing in g^ and g^ . 

The simplest renormalization factors [2]-[Z] depend solely on local hole 
densities x,- = 1 — n,-, i.e. 



9ij = \ f iiv^n 3tj= . ..w — -TT- 12 

^ y (xj + + 1) (xj + + 1) 

In the framework of our method, in the homogeneous case x = Xi, the addi- 
tional Lagrange multiplier coupled to the total particle number is introduced. 
For a non-homogeneous case, local chemical potentials should be introduced, 
which may be found in Refs. [5l El ITS] . 

In the case of more complicated form of g^, g^ factors, depending apart 
form Xi also on the inter-site mean-field variables ([31 El [III 112]), the MF 
Hamiltonian has to be modified further. The Lagrange multipliers related 
to the average hopping (bond order) and superconducting order parameters 
acquire non-zero values in the equilibrium situation. 

As pointed out previously, the form of Gutzwiller approximation of Ref. 
[3] does not reduce to multiphcation by renormalization factors, E.g., 
for the hopping amplitude, we have 



(Eqn. (15) of Ref. [3], but with different notation). This, in general, does not 
allow for identification of renormalization factor g^. However, for complete- 
ness, we want to compare the variational approach with the non- variational 
treatment of RMFT based on Ref. [3]. Hence, using fll3p we may quite rea- 
sonably define those factors for a simple homogeneous non-magnetic states 



^This corresponds to (0)0 in notation of Ref. [3] 
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[rii^ = riii = ni/2, (cJ^Cj^) = Xijt = Xiji = Xij) 



as 





— Ui 




-^i ( 




-f\ 




2 ^ 



(l-f)(l-!it)' 

The prescription f fTSj) will be referred to as a renormalization scheme (I). 
This will be confronted with (^-factors taken from Ref. [8] 



iTT\ „t ^ / AxiXj{l Xi){l Xj) , , 



4(l-x^)(l-Xj) 

(l-x2)(l-x2) + 8x.x,(|A.,f -|x.,f) + 16(|A,,r + |x.,r)' ^ ^ 

referred to as renormalization scheme (II). Usually within the RMFT, the 
physical (renormalized) superconducting parameter is defined as g^A instead 
of bare A itself, [2]. This is also the case for scheme (II). On the other hand, 
within the scheme (I), there appear a separate expression for renormalized 
value of (cj_o-Cjo-)G = {cj_aCia)Gi cf. Eqn. (18) of [3]. Again, for simple 
homogeneous, non-magnetic states, the corresponding factor can be quite 
reasonably identified as 



1'^ 




/;: 


( 




f\ 




2 ^ 



(I) 4 = Jf^Jf^(l + M;^7^)- (18) 



Note, that ([HD differs from ([H]). 



3 Results: d-wave superconducting and stag- 
gered flux phase 

Below we present our numerical results, first for the plain d-wave supercon- 
ducting state with no magnetic order (dRVB), and next for the staggered-flux 
(SF) normal solution. Even for those simplest MP states, the results obtained 
within various renormalization schemes and/or the solving methods differ re- 
markably. 
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3.1 d-wave superconducting (dRVB) solution at x = | 

The solution analyzed her^ is constructed to possess full symmetry of the un- 
derlying square lattice, only the superconducting order parameter is assumed 
to have dx2-y2 symmetry. Consequently, we are left with three independent 
mean fields A = (n, x. A), where Xx = X = Xy^ = ^ = ~^yi and the 
same number of the corresponding Lagrange multipliers, A = (A,A^,A'^), 
where A^ = AJ = A^, \^ = = -Xf, and n = rii = Xij = Xr, 

T = x{y) for the bonds in x (y) directions, respectively. Also, \/2Ajj = A^, 
and all the above quantities are taken as real. Diagonalization of in the 
present case yields 

Hx-fiN = J2 ^k(7lo7ko + + 5^ (a - E^) + C, (19) 



with 



= + (20) 
Dk = V2^D^cos{K), = -2^Tr cos{kr) - 1^ - X. (21) 

T T 

Also, 

Tr = -hrgl + \jr9iXr + A^, = \jr9i^r + Af , (22) 

X = + E {\jr9i{'2xl + A^) + 4x.A^ + 2A.Af ) . (23) 

T 

The Gutzwiller renormalization factors read now, respectively 



2-n \ (2-n)2 ; ' ^ {2 - 

A. . ^ 2(1 -n) / + 2A2\ , ^ 

,f(„.X,.A,) = ^(l + J^). (24) 



tTT\ tf N 2n(l-n) 

(//) gl{n,Xr)- 



2 ' 



n(2 -n)+ Axl 

(n, Xr, A,) = ^2(2 _ n)2 + (1 - n)2(4A2 - Sx^) + 4A4 + 16x^ ' ^^^^ 



*Part of the results presented in this section may be found also in |13) . 
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Note, that for (II) = g^. The generalized Landau functional is given by 

J-(1,A) = C(1,A) + J]((ek-i?k)-|ln(l + e-^^'<)), (26) 

k ^ 

from which the explicit form of Eqs. ([H]) can be easily obtained, but are not 
presented here in an explicit form . 

Numerical results. We begin with the analysis of dRVB solution for the 
'magic doping', x = | = 0.125. The equations ([8]) and ([9]) are solved for the 
lattice of A = A^Ay sites, A^. = = 256, with = 1 t^- = —3 J and for low 
temperature, /3 = 500. Both dRVB solution, as well as isotropic normal state 
(Fermi sea, (FS), not discussed explicitly here) are present. The staggered 
flux state, expected for lower doping, have not been found (var) or is not 
stable against FS (s-c) for x = 0.125. 

In Tables I. and II. we give the values of the thermodynamic potentials, 
mean-fields and molecular fields for dRVB solutions obtained within both 
methods {var, s-c) and both renormalization schemes ((I), (II)). 



Table I. Values of the thermodynamic potentials (per site) for dRVB 
solutions. Cl (F) stands for Q — XN {Q + fiN) for var and Qg-c 
{Qs-c + fJ's-cN) for s-c methods, respectively. 



Therm. Pot. 


var (I) 


var (II) 


s-c (I) 


s-c (II) 


Q/A 


-6.0444393 


-5.7586779 






n/A 


-1.0897421 


-1.0766359 


-1.01117344 


-1.03614582 


F/A 


-1.343195431 


-1.366146003 


-1.339864247 


-1.364716747 



For a fixed particle concentration n, the relevant thermodynamical po- 
tential is the free energy F. By construction of the solution the value of F 
obtained within each RS is alway lower for the var method then for the s-c 
one. The free energy is also the quantity, that determines which solution of 
([HD or OH]) corresponds to the stable equilibrium situation. However, by no 
means it may be used to favor one or another renormalization scheme. For 
example, if we compare the values of the F for var method in the present 
case, and also invoke its value F^^^., corresponding to the simplest Gutzwiller 
factors (given by Eqn. f|T2|) . and not analyzed explicitly here, c.f. however 
Ref. [12]), we see that = -1.5070 < F^^. < F^^^.. Clearly, it does not 
mean that RS defined by f[T^ should be preferred over (I) or (II). 

Table II. Values of the equilibrium chemical potentials and MF parameters 
{Aq, Ao) for dRVB solutions, /t stands for A + /i (var), and for ^s-c (s-c). 
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Variable var (I) var (II) s-c (I) s-c (II) 



-0.37551 
0.19074 

0.12345 

2.04556 
0.57209 

Let us now analyze Table II. First, within RS (II), |A^| is significantly 
larger than | A:^ | . This indicates that the Hamiltonian (fT9|) deviates from the 
Hartree-Fock form more with respect to Xt then with respect to A^. This 
is due to the particular form of Xr-dependent, A^-independent renormalized 
hopping term. Similar conclusions are valid for RS (I), but then the Xt- and 
At- dependence of (H) is more symmetric (c.f. Eqs. [TH and [T5|, or \2M and 
thus the difference between \X^\ and \X^\ is smaller then for RS (II). 




'r s M X s Y r 



Figure 1: (Color online) Dispersion relations along the main symmetry lines 
in the Brillouin zone for the dRVB solutions for a square lattice, of the size 
Ax = Ay = 256, and for the filling n = 0.875. Triangles - self-consistent, 
non- variational results for (I) and (II), red squares (I) and green circles (II)- 
the present variational method. 

Large values of |A^| significantly affect (through Tr and Dr, Eqn. fl22l) ) the 



fi 5.37285 

A -5.66251 

fi -0.28966 

Xx 0.19414 

A^ -0.15883 

^A, 0.10897 

^A^ -0.08915 

2T, 1.72718 

V2Dx 0.33830 



5.02004 
-5.35091 

-0.33087 -0.37565 

0.18807 0.19011 
-0.16985 

0.13200 0.12565 
-0.01111 

1.69732 2.01576 

0.58926 0.59565 
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quasi-particle spectra, E]^, cf. Fig[TJ Namely, the excitation energies within 
our method are always lower then those of BdG self-consistent approach. 
This is due to relative minus sign between equilibrium values of and Xr 
{X^ and At-), cf. also Ref. [13]. Although for RS (II) the differences between 
the methods {var, s-c) are pronounced mainly in the regions of the Brillouin 
zone which are far from the Fermi surface, the difference of the tangent at 
the cone near point S may be of some significance. However, what is more 
important, RS (I) var gives different excitation energies then other cases also 
along the X-Y direction, i.e. close to the Fermi surface (^k = 0). The reason 
for such behavior is obvious from the analysis of Tab. II; along X-Y direction 
■Ck = const and the main contribution to E]^ comes from Di^, determined in 
turn by D^, which is exceptionally low for (I) var. 

3.2 Doping dependence of the mean-field quantities: 
critical hole concentrations 




Figure 2: (Color online). Renormalization scheme (I): Doping dependence of 
the bond-order parameters Xx = Xyy the superconducting order parameters 
Aj, = — Aj^, and their renormalized correspondents g^Xx = 9^Xy ^-nd g^^x = 
—g^Ay, both for the s-c (triangles) and the var (squares) methods. 

Next, we are going to discuss the changes appearing as the function of 
doping. Those are the most interesting results obtained in the present paper. 

The (1 — n)-dependences of mean-fields Xt and A^^, as well as the phys- 
ical (renormalized) gap parameter g^A^ and renormalized hopping g^Xr for 
scheme (I) are analyzed in Fig. [2J (The analogous picture for RS (II) is 
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0.015 
0.01 
0.005 

°0 0.05 0.1 0.15 0.2 0.25 03 0735 
X = 1-n 

Figure 3: (Color online). Doping dependence of renormalized superconduct- 
ing order parameter g^/S.^ for the paired (dRVB) state within various schemes 
described in the text. For x = 0.2, from the top to the bottom: (I) s-c (blue 
triangles), (II) var (green circles), (II) s-c (violet triangles), (I) var (red 
squares) . Only the latter choice yields the (right value of the) critical doping 
Xc ~ 0.27, above which the gap disappears. 

given in f3]). Note, that for (I) var, vanishes at the critical concen- 
tration Xc ~ 0.27. This is in much better agreement with the experimental 
results then the predictions of the other cases, for which Xc > 0.35. The 
situation is illustrated explicitly in Fig. El where the x-dependence of the 
g^Ax is shown also for the scheme (II), for both methods. Also, the right 
value of the optimal hole concentration {x ~ 0.125) is obtained for (I) var, 
in contrast to either non- variational treatment or to RS (II). 

3.3 Staggered flux solution at x « | 

As a next example we analyze the staggered-flux (SF) phase. This MF 
state has a long history, first being proposed by Affleck and Marston [22], 
as a variational trial MF state for the Heisenberg model. It was intensively 
investigated later, due to its possible connection with the pseudo-gap state 
in cuprates, cf. e.g. [23] . 

The SF differs from a normal Fermi sea (FS) solution by the presence 
of complex hopping amplitude Xij = \x\ exp({—lY^'^~^^'^Hip) = ± i^2- Such 
Xij implies existence of circulating currents, which direction changes from 
plaquette to plaquette in an alternating fashion (orbital antiferromagnet), 
[23] • Consequently, a two-sublattice structure emerges, with the unit cell of 
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the size a in direct space and new (folded) Brillouin zone (NBZ). 

Within the framework of our method we add appropriate constraints. 
This introduces, apart from A ascribed to n, also the complex Lagrange 
multiplier rjij = rji^ iri2, tailored to Xij (our sign convention for rjij for each 
bond is opposite to that for Xij)- Thus, we have three independent real 
mean fields A = {n, ^1,^,2), and the same number of the corresponding real 
Lagrange multipliers, A = {X,i]i,ri2). In the present case the renormalization 
factors (fT5l) and (fTTI) read, respectively 

/TT^ f, ^ 2n(l — n) 

{11) 9\n,^,,^2'- 



n(2-n)+4(ef + ei)' 

^ = n^i2 - nf - 8(1 - nfUl + H) + ml + ^1)^ ' ^^^^ 

Instead of Eqs. f fT9|) - ([7]) we have now 

H^-^N = C{A, X)+Y, ^w«L«ks., (29) 

ksCT 



where Eks^ = -/i + s^/e^ + xl, Ck = Ti7+(k), and Xk = T27_(k). Also, 
7-1- (k) = 2(cos(A;xO ± cos(A;j^)) and 

Ti = tg' + pg'^^ +771, T2 = ^ J^?-'^ + V2, (30) 
C = A (An + 3(7^(^2 _^ ^2) ^ 8(^^^^ ^ ^^^^^^ _ ^3^) 
The generalized Landau functional ([7]) takes now the following form 

^ NBZ 

J^iAX) = C(AA)--5^1n(l + e-'^^-^). (32) 

Numerical results. As mentioned in Sec. ??, at x = 1/8 the existence of 
SF solution of Eqs. (|8]) has not been numerically confirmed (for var method), 
whereas for the s-c method the SF solutions of of Eqn. (Q have been found 
unstable against Fermi sea (FS) with ^2 = 0. In all four cases the SF— J-FS 
transition is located at the critical concentrations 0.11 < Xc < 0.12. How- 
ever, our numerical procedures for var method turned out to be unstable 
in the vicinity of the Xc- For that reason, we chose doping x = 13/128 
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{n ^ 0.898), which is in a safe distance from each of Xc, but for which the 
differences between the methods are pronounced (as they generally increase 
with increasing doping). 

The parameters of the Hamiltonian are the same as for the RVB case, 
except sign convention for t, now t = 3 (t = — 3 in the RVB case), also 
Arc = Ay = 512. Again, we work with low T = 1/500. The thermodynamic 
potentials and mean-field variables are listed in Tabs. Ill and IV. 

Table III. Equilibrium values of the thermodynamic potentials (per site) for 
SF solutions, for n ^ 0.898, = 512. (F) stands for n - XN {Vi + /dN) 
for var and flg-c {^s-c + fJ's-cN) for s-c methods, respectively. 



Therm. Pot. 


var (I) 


var (II) 


s-c (I) 


s-c (II) 


n/A 


-5.90103948 


-5.75555848 






n/A 


-0.69003640 


-0.71583672 


-0.49006797 


-0.49344983 


F/A 


-1.18762800 


-1.16898668 


-1.18536044 


-1.16599431 



Table IV. Values of chemical potentials and MF parameters for SF solutions 
for n ^ 0.898, A^ = 512. fl stands for A + (war), and for /i^.^ (s-c). 



Variable 


var (I) 


var (II) 


s-c (I) 


s-c (II) 




5.24623 


5.10505 






A 


-5.80007 


-5.60943 






P- 


-0.55384 


-0.50438 


-0.77389 


-0.74857 


ei 


0.19222 


0.19321 


0.18805 


0.18844 


6 


0.10298 


0.09840 


0.11856 


0.11731 


Vi 


-0.13476 


-0.14848 








-0.07220 


-0.07562 






Saf 


-0.23514 


-0.22522 


-0.24436 


-0.23527 


$n 


0.31311 


0.29987 


0.35811 


0.35449 


Ti 


0.80695 


0.77919 


0.92799 


0.91123 


T2 


0.18241 


0.16008 


0.29312 


0.28009 



Within the standard mean-field approach (e.g. [12]) the (fictitious) flux is 
defined as $□ = 2^ X](jj)Gn ^^g(Xji)) and □ denotes plaquette composed from 
four bonds. Also, for s-c method, i.e. for rji = ri2 = 0, we have 



Arg(x 



arctan | — 



arctan 



To 



Ti - tg^ 



(33) 



the last equality follows from f l30|) . Interestingly, this equality holds also for 
the variational approach with rji 7^ 0, 772 7^ 0. It can be shown analytically, 
that ^1/1^2 = V1/V2, which together with (I5U|1 yields Thus the two 
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possible and a priori different definitions of $□ within var method turn out 
to be equivalent. Antiferromagnetic correlations are defined on the MF level 

as saf = -I9' {ei+e2), mm- 




Figure 4: (Color online) Dispersion relations for both upper (-Ek+) and lower 
(-Ek-) subbands for the SF solutions along the main symmetry lines in the 
Brillouin zone of the square lattice, of the size = A^^ = 512, and for the 
filling n ^ 0.898. Triangles - self-consistent, non-variational results for (I) 
and (II), red squares (I) and green circles (II)- the present variational method. 
Explicitly, for (-E'k+) near the maximum at point M, from the bottom to the 
top: (II) var, (I) var, (II) s-c (violet triangles), (I) s-c (dark blue triangles). 

In contrast to the dRVB case, now the differences between renormaliza- 
tion schemes ((I) vs (II)) within each {var, s-c) method are small. This 
is because in the absence of Ajj, the Xij -dependences of g^^j-^ and gj^jjs^ are 
quite similar, and the Xij -dependence of g^jj-^ is weak for doping x ~ 0.1, 
thus causing no quahtative and only minor quantitative differences between 
the renormalization schemes. On the other hand, the generic modifications 
introduced by the variational approach within each RS are more significant. 

Namely, from Tab. IV we see, that s-c method favors SF more strongly 
than var method, which is indicated by the values of $□ and ^2- Also, Ti 
and T2, determining quasi-particle spectra, are smaller within var method, 
and so are the quasi-particle energies, c.f. Fig. HI 
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4 Concluding remarks 



In summary, in this work we have compared, within two methods of approach, 
the two Gutzwiller renormahzation schemes for the renormahzed mean-field 
theory (RMFT) of t-J model in its simplest form. We emphasize the advan- 
tages of the combination of the renormahzation scheme of Ref. |3j with the 
variational method proposed in Ref. [2]. First, a number of theoretical ar- 
guments strongly favor this choice. Moreover, in contrast to the other cases 
investigated by us (e.g. either non-variational method or renormahzation 
scheme of Ref. ^), the former approach quite correctly predicts the upper 
critical doping Xc ~ 0.27 for a disappearance of the SC order. Also, the value 
of the optimal doping x ~ 0.125 is quite correctly predicted. In our opin- 
ion, the formalism of Ref. p] augmented with the self-consistent variational 
treatment, gives a chance for the complete and consistent one-particle de- 
scription (in the form of RMFT) for a t-J model. Such description, however, 
must encompass the t-J model in its complete form [I], and include also more 
complicated symmetry breaking patterns. This is the subject of our current 
investigation. 

Appendix A: Equivalence relation for mean-field Hamil- 
tonians 

Below we present some details of our formalism, which are necessary for the 
present dissension. We also comment on the relationship between our method 
and the formalisms of Ref. [10] . We start from the MF Hamiltonian of the 
form 

M M' 

H{A) = Heo + Co{A) -Id + YI Cs{A)A, + G^(l)5^. (34) 

s=l w=l 

In the above, Heo is an A-independent part, Co{A), Cs{A) and G^iA) are 
some complex-valued functions of mean-fields A. Operators are those, 
which average values are not present in ^^(^4). We also assume that all 
the operators appearing above are bilinear in creation and/or annihilation 
operators. From flMl) we have 

M 

Hx{A) = H{A)-J2U^s-As) 

s=l 

M M' 

= H,o + C^{A) -Id + YI C^i^)^s + Y GUA)B^, (35) 

s=l w=l 
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with Cq = Co + J2sLi '^s^s, Cg = Cs — As. For a given value of A, the self- 
consistency equations (second half of the 2M equations ([H])) niay be written 
as 



At = {At)x = TT[Atpx] 

Tr[A exp ( - /3(Keo + EsU C^iMs + T^Li G^{A)B^))] 



Tr[exp ( - /3(i^eo + Y.ti CKA)As + Ei=i G^{A)B^)] 



(36) 



In above, the Cq{A) term canceled, and K^q = Heo — l^N . Eqs. fl5Bl) may 
be formally solved (which usually cannot be achieved in an analytic fashion) 
for A, and then A = A (A), thus for a given A, Eqs. fl36|l determin^ Cg{A). 
As a result, the A-dependence of the coefficients C'^{A) is determined by the 
choice of operators ^feo, and {1?^,}^^^ If we change the original 

form of the coupling of {As}^^i operators to the mean fields, according to 

C,{A)-^Cs{A)=^s{C,{A)), (37) 

but without changing Gu){A) functions, then the Lagrange multipliers change 
in a way that Cg{A) is unchanged, in order to fulfill (l36l) . The only A- 
dependent part which may be non-trivially modified by (1371) is Cq{A) ■ Id- 
Obviously, S\, defined as (c.f. [H]) 



r . . . 1 1 

5a = Tr -px\nf)x-l3{p\H-p.p\N-^^Xsh{^s-As))-u{px-—) (3 



D' 

3 = 1 

and J-'{A,\) = —f3~^S\{px{A,\),A,\) may be also modified. However, the 
density operator px is invariant under (157|) . so are, for given A, all the aver- 
ages, also those of Heo and Bu, operators. Moreover, if the transformations 
fl37|l are such that {Hx)x = {H)x remains unchanged, the value of J^{A, X{A)) 
is not modified. Consequently, in such a situation, the equilibrium values of 
mean fields, as well as of the quantities C^{A) (e.g. T^-.D^ of Eqn. fl22]) 
or Ti,T2 of Eqn. flHT]) ) are also invariants of ([17]). Summarizing, all MF 
Hamiltonians constructed from the same set of operators Heo, {As}gL^ and 
{Bj^,}^^^ and having the same A-dependence of the expectation value, are 
equivalent, and belong to the same equivalence class. The transformations 
fl37|) may be viewed in analogy to gauge transformations, not changing the 
physical content of the model. 



^We assume that the solution exists and is unique, or there exist a finite number of the 
solutions, but one of them may be unambiguously selected. 
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Suppose, that the MF Hamiltonian flM]) is such that 

i^eo = 0, V^: G^ = 0. (39) 

Then the equihbrium values of may be easily obtained in an analytic 
fashion, using Eqn. (17) of Ref. [13], i.e. 

Now, we choose a transformation (!37|) of a specific form, 

Co{A) ^ (7o(^) = Co(l) + ^ C,(l)^. (41) 



This yields 



M 



Hi^\A) = -J2 - A,) + W{A), (42) 



s=l 



with H{A) = W{A)iD, W{A) = {H{A)) = Co{A). Using (gO]) we may 
rewrite (H2l) as 



^r^(^) = E ^^(^^ - + (43) 

s=l 



because avr(A)/M, = d{H{A))/dAs = {dH{A)/dA,). The form (USD is the 
most convenient, as half of the variables (A) are eliminated, which reduces 
the number of equations to be solved numerically. 

Also, in the limit {(3 — )■ oo) the results of finite-temperature formalism are 
essentially identical to those the true T = analysis, and we may compare 
them with those of Ref. [TU]. Hamiltonian of that Reference reads in our 
notation 

Hi'<\A) = -j:X.(A)A. = j:^A., (44) 

s=l s=l * 

and differs from (l42l) only by the constant term X^fli ^sAs + {H{A)). Conse- 
quently, both ( 142|) and fH4|) have the same eigenvalues and eigenvectors (the 
presence of —fiN obviously does not change the above arguments). Please 
note, that the analytical evaluation of Lagrange multipliers through (HUIl . 
and hence the application of the method of Ref. [TU] is possible only if the 
conditions fl39|l are fulfilled, but not for the general form flMl) of the MF 
Hamiltonian. 
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